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ABSTRACT 


A  method  is  developed  for  calculating  the  mass  flux 
to  a  liquid  droplet  surrounded  by  its  pure  vapor  as  a 
function  of  the  Knudsen  number.  The  Knudsen  number  K  is 
defined  as  the  ratio  of  the  mean  free  path  to  the  droplet 
size.  When  the  mass  flux  is  expanded  in  terms  of  the 
inverse  Knudsen  number  a  -  K-1,  we  obtain  a  series  of 
the  form 


r  =  r 


(0) 


r(1V  + 


f(2)a2lna 


It  is  shown  that  the  coefficients  are  determined  by  inte¬ 
grals  associated  with  sequences  of  successive  collisions 
among  a  number  of  vapor  molecules  and  the  droplet.  In 
particular,  we  derive  the  collision  integrals  for  the 
first  three  coefficients  of  the  inverse  Knudsen  number 
expansion  for  r.  These  collision  integrals  bear  a  close 
similarity  to  the  collision  integrals  derived  in  earlier 
technical  reports  for  the  density  dependence  of  the 
transport  properties  of  gases.  It  will  be  demonstrated 
in  a  subsequent  technical  report  that  the  same  method  can 
be  used  to  calculate  the  aerodynamic  force  on  an  object 
in  a  gas  stream  as  a  function  of  the  Knudsen  number. 
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For  a  droplet  to  be  in  equilibrium  with  the  surrounding 
vapor,  its  rate  of  vaporization  must  be  equal  to  the  rate 
at  which  vapor  is  condensing  onto  the  surface  of  the 
droplet.  Such  an  equilibrium  is  only  possible  when  the 
droplet  has  a  very  particular  critical  size  [1] .  Droplets 
smaller  than  the  critical  size  are  called  "embryos1';  they 
will  tend  to  evaporate  and  disappear.  Droplets  larger 
than  the  critical  size  will  tend  to  grow.  Nucleation 
theories  are  concerned  with  predicting  the  rate  at  which 
these  growing  nuclei  are  created  per  unit  volume.  For  a 
review  of  the  theory  of  nucleation  the  reader  is  referred 
to  reference  [2]  . 

In  this  technical  report  we  are  concerned  with  the 
rate  of  change  of  the  liquid  droplets  once  formed  in  the 
vapor.  In  particular  we  shall  consider  the  rate  of  change 
of  mass  of  liquid  droplets  surrounded  by  the  pure  vapor. 

The  theory  to  be  presented  is  valid  for  droplets  formed 
either  by  homogeneous  or  heterogeneous  nucleation. 

We  assume  that  the  average  distance  between  all  liquid 
droplets  is  significantly  larger  than  both  the  radius  R 
of  the  droplets  under  consideration  and  the  mean  free  path 
Z  of  the  molecules  in  the  vapor  phase.  Under  these  circum¬ 
stances  we  may  neglect  the  interaction  between  different 
droplets  and  it  is  sufficient  for  the  problem  of  droplet 
growth  to  evaluate  the  net  mass  flux  of  vapor  molecules  to 
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an  individual  droplet  of  radius  R  and  temperature  Tx , 
surrounded  by  its  vapor  at  temperature  T  .  When  the  density 
of  the  vapor  is  not  too  large,  this  net  mass  flux  can  be 
obtained  by  solving  the  Boltzmann  equation  in  the  vapor 
phase,  subject  to  the  appropriate  boundary  conditions  at 
the  surface  of  the  droplet. 

The  mass  flux  to  the  droplet  will  in  general  depend  on 
the  ratio  a  =  -j  (the  inverse  Knudsen  number)  .  In  the  free 
molecular  limit,  a  -*■  0,  this  mass  flux  is  given  by  the 
well  known  Hertz-Knudsen  formula  which  is  reproduced  in  the 
next  section.  For  a  liquid  droplet  surrounded  by  a  mixture 
of  its  vapor  and  an  inert  gas,  the  mass  flux  in  the  contin¬ 
uum  limit  is  given  by  Maxwell’s  equation  [3]  .  Only  a  few 
investigators  have  studied  the  problem  in  the  intermediate 
regime.  Brock  [4]  has  considered  the  first  correction  term 
to  the  free  molecular  flux  for  a  binary  gas  mixture,  but 
replaced  the  Boltzmann  equation  with  a  linearized  BGK 
model  equation.  Shankar  15]  and  also  Sampson  and  Springer 
[6]  approximate  the  distribution  function  by  a  two  stream 
Maxwellian;  the  adjustable  parameters  in  this  function  are 
then  determined  so  that  the  distribution  function  will 
satisfy  a  limited  number  of  moments  of  the  Boltzmann 
equation.  These  approaches  each  involve  approximations 
whose  nature  and  magnitude  are  not  clear  and  which  there¬ 
fore  make  it  difficult  to  estimate  the  range  of  their 
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validity.  It  would  thus  appear  to  be  useful  to  develop  a 
more  rigorous  solution  to  the  problem.  The  research 
reported  here  is  an  attempt  in  that  direction. 

We  shall  formulate  a  systematic  expansion  of  the  mass 
flux  to  a  liquid  droplet  surrounded  by  its  pure  vapor  in 
terms  of  the  inverse  Knudsen  number  a.  This  approach  will 
be  adequate  in  the  nearly  free  molecular  regime,  where  a 
is  small  compared  to  unity.  This  condition  is  encountered 
in  many  situations  involving  either  homogeneous  or  hetero¬ 
geneous  nucleation,  where  the  radius  of  the  critical  size 
clusters  is  indeed  considerably  smaller  than  the  mean  free 
path  [7],  When  the  rate  of  change  of  mass  of  the  droplet, 
r,  is  expanded  around  the  free  molecular  limit  (a->0)  ,  we 
obtain 


,<0) 


+  r(1)a 


+  f(2)ct2lna  + 


(1-1) 


The  successive  coefficients  of  this  series  will  be  given 
by  integrals  associated  with  sequences  of  successive  colli¬ 
sions  among  an  increasing  number  of  vapor  molecules  in  the 
presence  of  the  droplet.  We  shall  in  particular  derive 
explicitly  the  collision  integrals  for  the  first  three 
coefficients  in  the  expansion  (1-1). 

The  structure  of  the  expansion  (1-1)  is  analogous  to 
the  density  expansion  of  the  transport  properties  of  gases 
[8].  In  fact,  we  shall  show  that  the  collision  integrals 
for  the  mass  flux  have  a  close  similarity  with  the  collision 


4 


AEDC-TR-72-172 


integrals  derived  in  earlier  technical  reports  for  the 
transport  coefficients  of  a  moderately  dense  gas  [9,10]. 

A  similarity  of  the  same  kind  between  the  density  expansion 
of  the  transport  coefficients  and  the  Knudsen  number  expan¬ 
sion  for  the  drag  coefficient  of  objects  in  rarefied  gas 
flows  was  previously  noted  by  Dorfraan  et  al.  [11,12]. 

In  Section  II  we  shall  formulate  the  problem  and 
specify  the  boundary  conditions  employed  in  the  analysis. 

In  Section  III  we  shall  derive  a  set  of  well-defined  inte¬ 
grals  that  determine  the  coefficient  of  the  first 

correction  term.  In  the  same  section  we  develop  a  diagram¬ 
matic  notation  to  represent  sequences  of  collisions  between 
the  molecules.  This  notation  will  enable  us  to  elucidate 
how  the  integrals  of  the  expansion  (1-1)  for  r  are  related 
to  specific  sequences  of  collisions  among  the  molecules  in 
the  presence  of  the  droplet.  In  Section  IV  we  shall  show 

that  the  next  term  is  logarithmic  in  the  Knudsen  number 

~  (2) 

and  derive  collision  integrals  for  its  coefficient  r  . 

In  Section  V  we  discuss  the  relationship  between  the  new 
collision  integrals  and  those  derived  earlier  for  trans¬ 
port  coefficients  and  drag  coefficients. 
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II.  Steady  State  Mass  Flux  to  and  From  a  Liquid  Drop 

We  shall  calculate  the  net  mass  flux  to  a  spherical 
liquid  drop  with  radius  R  and  temperature  T^  surrounded  by 
its  pure  vapor.  The  coordinate  system  is  chosen  so  that 
the  center  of  the  droplet  is  located  at  the  origin.  The 
number  density  and  temperature  of  the  vapor  at  large  dis¬ 
tances  from  the  droplet  are  maintained  at  n^  and  T  , 
respectively. 

We  shall  follow  the  previous  authors  [4,5,6]  by  consid¬ 
ering  the  rate  of  change  of  the  number  of  molecules  in  the 
droplet  in  the  steady  state  approximation.  That  is,  we 
assume  that  the  droplet  is  sufficiently  large  so  that  we 
can  neglect  the  change  in  the  radius  and  temperature  of 
the  drop  when  a  limited  number  of  molecules  condense  onto 
or  evaporate  from  it.  Then  the  single-particle  distribution 
function  F(r,v^)  of  the  molecules  in  the  vapor  phase  will 
satisfy  the  time- independent  Boltzmann  equation 

=  fd^2fd^i2k12  [I'Cr,v1,)F(r,v2’)-F(r,v1)F(r,v2)  J, 

(2-1) 

where  r  is  a  vector  from  the  droplet  center  to  some  point  in 

A 

the  vapor,  d^2  the  perihelion  vector  specifying  the  geometry 
of  a  collision  between  two  molecules  and  k12 ^21*^12^  a 
positive  scalar  factor,  proportional  to  the  differential 
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cross  section  for  a  collision,  in  the  notation  of  Chapman 
and  Cowling  [13]  .  The  velocities  are  related  to  v^Vg 

by: 


■y 

V 


1 


•6 


12' 


C'l2)012» 


(2-2) 


-y  ->  -> 

where  v  .  .  =  v .  -  v . . 
i  J  i  J 

We  need  to  solve  the  Boltzmann  equation,  (2-1),  subject 
to  the  appropriate  boundary  conditions  both  at  the  surface 
of  the  droplet  and  at  infinity.  Far  from  the  droplet  the 
pure  vapor  will  be  in  equilibrium,  so  that: 


lim  F(r ,  v) 


(2-3) 


where  m  is  the  mass  of  a  vapor  molecule  and  k  is  Boltzmann's 
constant . 

We  have  already  assumed  that  the  liquid  droplet  is  a 
sphere  with  radius  R.  To  formulate  the  boundary  conditions 
at  the  droplet  surface,  we  must  specify  the  mechanism  for 
the  interaction  of  the  vapor  and  the  liquid  surface.  This 
interaction  is  usually  described  in  terms  of  a  mass  accom¬ 
modation  coefficient.  For  simplicity,  we  shall  follow 
Shankar  [5l  and  take  this  accommodation  coefficient  to  be 
unity,  i.e.  all  molecules  that  strike  the  surface  will 
stick  to  it.  All  molecules  that  leave  the  surface  shall 
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be  assumed  to  be  emitted  diffusively  with  a  Maxwellian 
distribution  having  the  droplet  temperature  and  a 
number  density  nL  corresponding  to  the  saturation  pressure 
P  =  n.kTT  about  a  droplet  of  radius  R.  Thus  the  boundary 
condition  on  the  distribution  function  at  the  droplet 
surface  is: 

F(R,  v)  =  nL  (j”  )2exp  (-  g*  )  for  ?.1>0.  <2-4) 

L  Li 

The  pressure  ,  and  thus  n^,  is  related  to  the  saturation 

vapor  pressure  P  for  the  vapor  in  equilibrium  above  a 

p«  s  * 

planar  surface  of  its  liquid  at  temperature  T  via  the 

Lt 

Thomson-Helmholtz  equation  [ 1]  . 

We  note  that  (2-4)  represents  the  distribution  function 

at  the  droplet  surface  of  molecules  in  the  vapor  phase. 

Thus,  strictly  speaking,  F(R,v)  is  to  be  interpreted  in  this 

report  as  lim  F(R+£,v).  With  this  in  mind,  the  net  mass 
6^-0 

flux  T  to  the  spherical  droplet  may  be  written  in  terms  of 
the  distribution  function  as: 

T  =  -47;R^m  ^dv  v*R  F(R,v)  ,  (2-5) 

where  R  =  S/R  is  the  unit  vector  in  the  direction  of  R.  The 
number  of  molecules  in  the  droplet  will  increase  or  decrease 
depending  on  whether  r  is  positive  or  negative.  Since  the 
net  mass  flux  obviously  vanishes  when  nT  =  n  and  Tt  =  T  , 

Lj  OD  L  CO7 
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it  becomes  our  task  to  evaluate  r  when  n.  j*  n  and/or 

J j  00 

Tt  /  T  . 

The  purpose  of  this  report  is  to  elucidate  the  struc¬ 
ture  of  the  collision  integrals  that  determine  the  variation 
of  r  with  the  Knudsen  number.  The  assumptions  that  the  drop¬ 
let  has  a  spherical  shape  and  that  the  mass  accommodation 
coefficient  is  unity  are  not  essential  to  our  analysis,  but 
are  introduced  for  simplicity.  The  collision  integrals  to 
be  derived  can  also  be  formulated  for  droplets  of  any  geo¬ 
metry  and  the  method  can  be  generalized  to  the  case  where 
the  mass  accommodation  coefficient  is  smaller  than  unity. 

In  order  to  exhibit  the  expansion  parameter  in  the 
perturbation  procedure,  it  is  convenient  to  introduce 
dimensionless  quantities  [14] .  For  this  purpose  we  define: 


P  R 


03 


/  2kT  A 
-  ' 


I  = 


CD. 


12 


* 

n  = 


k12 

2  ' 

7T0  03 

no 

n. 


n 


3 


03 


f  (p,t)  =  ~  F(r,v) 


*  L 
T  a  • 

T 

00 


(2-6) 


where  o  is  a  length  parameter  characterizing  the  range  of 
the  interaction  potential  between  the  molecules.  The  time- 
independent  Boltzmann  equation  may  then  be  written  in 
dimensionless  form: 
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?l  •  -  a/d?2/d8 12*12  [f(p.l 

(2-7) 

o 

with  a  =  7T0  Rn^.  (2-8) 

The  boundary  conditions  (2-3)  and  (2-4)  become: 

-*  "*■  —  2-  n 

lira  f  (p,£)  -  7T  *  exp(-(-  )  for  all  (2-9a) 

p-fco 

f(p,£)  =■  n  (ttT  fTexp(-  for  £-p>0  .  (2-9b) 

T 

and  the  net  mass  flux  r  to  the  droplet  is  given  by: 

T  =  ~47TR2mnooaioo  jfd£  £*p  f(£,|)  .  (2-10) 

Since  the  mean  free  path  i  of  the  vapor  molecules  is 
of  the  order  1  ArC  nw,  the  dimensionless  parameter  a  in  (2-7) 
is  indeed  the  inverse  Knudsen  number  R/J.  In  this  report 
we  consider  the  nearly  free  molecular  regime  a<<l  and 
solve  the  Boltzmann  equation  by  a  perturbation  procedure  in 
which  we  consider,  successively,  terms  of  higher  order  in  a. 
Unfortunately,  the  name  "Knudsen  number  iteration”  is  some¬ 
times  associated  with  a  procedure  in  which  one  attempts 
to  represent  the  solution  of  the  Boltzmann  equation  by  a 
power  series  in  a  [ 15] .  Such  a  procedure  leads  to  diver¬ 
gence  difficulties  which  have  the  same  physical  origin  as 
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the  divergences  encountered  when  the  transport  coefficients 
of  a  gas  are  represented  by  a  power  series  in  the  density 
l8, 16].  As  already  pointed  out  by  Willis  [17],  it  is  not 
necessary  in  the  perturbation  procedure  to  impose  the 
assumption  that  the  solution  be  a  power  series  in  a.  We 
shall  follow  a  modified  Knudsen  number  iteration  procedure 
in  which  we  collect  terms  that  are  of  successively  higher 
order  in  a,  but  do  not  prejudice  the  result  by  assuming 
that  these  terms  are  all  powers  of  a. 

The  free  molecular  solution  f^(p,f)  is  obtained  in 
the  limit  a  ■+  0,  so  that: 

t  •  ~  )  -  o  .  (2-n: 


This  equation,  together  with  the  boundary  conditions  (2-9), 
implies  that  f^(p,£)  is  a  two-stream  Maxwellian  given  by: 


f 


(0) 


7T  2  exp(-£2) 
n*(7rT*)  rexp(- 


for  |  e  I(p) 
for  £  c  II  (p) 


(2-12) 


For  any  p,  the  regions  I (p)  and  II (p)  are  velocity  regions 
defined  by: 


Region  I 


-+•  •* 

^P'%<  0  or 

tp.?>0  and  |px||  2  -  e2  >  0  , 


(2- 13a) 
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Region  II  p-|>0  and  |pxf|2  -  %2  <  0  .  (2-13b) 

For  a  given  p,  the  zeroth  order  term  f  (p,£)  is  clearly 
discontinuous  in  velocity  space.  Region  II  is  the  free 
molecular  cone  of  influence  associated  with  the  droplet  as 
illustrated  in  Fig.  1;  molecules  having  a  velocity  in  this 
region  have  originated  from  the  droplet  and  have  arrived 
at  the  position  p  without  suffering  any  collision.  Mole¬ 
cules  having  a  velocity  in  region  I  are  unaffected  by  the 
presence  of  the  droplet  and  may  be  considered  as  having 
arrived  from  infinity  by  free  streaming. 

As  a  short  hand  notation  we  shall  write  the  two-stream 
Maxwellian  (2-12)  as: 

f(0)  <P,f>  =  f£0)<p,f)  +  fj [^(p,?)  ,  (2-14) 


where 


fTco)Cp,t)  = 


7T 


-  2 


exp  (-£2)  -  fm(£)  for  i  e  I(p) 


0  for  \  e  1 1  (p)  , 


(2-14a) 


and  = 


[n*  (7rT*)  2  exp(-  ^-*)  for  %  e  II  (p) 

■t  T 

0  for  |  e  I  (p)  .  (2- 14b) 


12 


Fig.  X.  The  velocity  regions  I  (p)  and  II (p)  at  a  given  position  p  in  the 
vapor.  Region  II  (p)  is  the  free  molecular  cone  of  influence 
associated  with  the  droplet. 
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The  net  mass  flux  r  to  the  droplet  in  the  free  molecular 
approximation  reduces  to: 


which  is  the  familiar  Hertz-Knudsen  formula  [181. 


-15) 
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III.  First  Correction  Term  to  the  Free  Molecular 
Approximation 

m  -*• 

In  order  to  obtain  the  first  correction  terra  fv  '(p,|) 
to  the  free  molecular  distribution  function  we  approximate 
f(p,f)  in  the  right  hand  side  of  the  Boltzmann  equation 
(2-7)  by  the  free  molecular  solution  f^(p,£).  We  thus 
consider: 

x  .sf(1)<p.ii> . 

5l  £ 

-  aJdI2ydS12r12[f (0>  (p,  5 p  t (0)  (0)  0, % 1 (0>  <jS,f2>] 

(3-1) 

Since  the  free  molecular  solution  (2-12)  satisfies  the  boun¬ 
dary  conditions  (2-9)  exactly,  we  require  that  all  higher 
order  approximations  f ^ (p,f)  for  i>  0  satisfy  the  boundary 
conditions: 


(i)(p,|)  -  0 

for  all  |  , 

(3-2a) 

U)  Cp,b  -  o 

for  l*-p>0  . 

(3-2b) 

Equation  (3-1)  can  be  integrated  to  yield  [14]  : 

t 

f  (£>?!>  -  f  ^  (P-Ijt,^)  +  ,1^  ,  (3-3) 

0 


15 


AEDC-TR-72-1  72 


^  T 

where  J(p,£)  represents  the  right  hand  side  of  (3-1).  In 
order  to  make  use  of  the  boundary  conditions  (3-2) ,  we 
select  the  limit  of  the  integration  such  that: 


t 


e  Up) 

for  e  II  (p) 


(3-4) 


where 


/(p-ij)2  + 1  -  p2 


(3-5) 


is  the  time  it  would  take  a  molecule  with  velocity  ^  e  II  (p) 
to  travel  from  the  droplet's  surface  to  the  position  p. 

This  results  in  the  vanishing  of  the  integration  constant 
in  (3-3)  for  the  problem  under  consideration.  We  thus 
obtain  for  the  first  correction  term: 


f(1)(p>I1)  -  +  f ^(P.tj) 


(3-6a) 


with 


fI  )  (^’^1)  afdT12fd^2j^l2K12[f  C0)  (p“^lT12’^l)f  (0)  (^"^1t12,^2) 

0 

-  f  (0)  (P_?1T12’^i)f  for  ^1  €  I(p) 


f 


(l) 

I 


-  0 


for  e  II(p)  i 


(3-6b) 
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and 


T(p,ei) 


^/'dri^/d^2Jd®12K:12  5  ^“^lT12'^l)f  (  ^P*"^  iTi2  2^  + 


-  f(°)<P~1iTi2>?i>f(0)^-ilTl2^2)l  for^l  £  11  <P>* 


fII)  <P’?l>  "  0  f°r  h  C  I(P> 

<3-6c) 

The  T12“in^eSration  in  these  integrals  extends  over  the 
free  trajectory  traversed  by  a  molecule  1  with  velocity 
If  we  substitute  the  expansion  f(p,£)  =  f^(p,^)  + 

+  f^(p,£)  +  .  .  .  into  the  expression  (2-10)  for  the  mass 
flux  to  the  droplet,  we  obtain: 


r  =  r(0)  +  a  r(1)  +  .  .  .  , 


(3-7) 


with 

<D  _ 


=  -47rR2mnoo(jjooa"1  jd£j  fj/P  f (1)  (p.lj) 


“  -^mn^  jTdr  12  /d312K12  \rp 

0 

4 1  ■  p<  o 

f(0)<P-«ri2^i)fC0)(P-?lT12'«2V 


-1  (0>  Cp-iXT12,t{>  f  (0)  (p-?lTl2,I2> 


(3-8) 
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The  first  correction  term  to  is  proportional  to  the 

inverse  Knudsen  number  a  provided,  of  course,  that  the 
integral  (3-8)  is  convergent. 

In  this  report,  we  assume  that  the  interaction  poten¬ 
tial  between  a  pair  of  vapor  molecules  has  a  finite  range 
which  is  much  smaller  than  R.  It  is  then  possible  to 
decompose  into  a  set  of  convergent  integrals  each  of 

which  is  related  to  a  particular  sequence  of  successive 
collisions.  For  this  purpose  we  use  the  relation 

f^  (£,  i')  f^  =  f00^i>foo^2)  *or  the  distribut  ion  function 

at  infinity  and  rewrite  r  ^  as: 

00 

r  =  -47TR  mn0Oa5o^/dT i gjd £  ^  1 2*'  1 2  ^l’p 

°  ?T-P<0 

.  [  f  (0)  0-1  ir12,X{>  f  (0)  (P-ifw  ?a') -f„  (4  {)  f M  C€2o]  + 
00 

-  4-nR2mDjaJaT1Jat1iiJ'ai12^12  ?rp 
°  Vp<° 

•[f„«l>^(«25-*(0)  <P-1lT12' <0>  0-?iT12.t2)]  ■ 


(3-9) 

As  a  next  step  we  separate  this  integral  into  terms  asso¬ 
ciated  with  the  contributions  f^0-*  and  f^  of  the  free 
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molecular  flow  solution  f^.  For  this  purpose  it  is 
convenient  to  distinguish  in  f  contributions  from  inside 
and  outside  the  free  molecular  cone  of  influence.  We 
therefore  write: 


f^te)  =  f“(p,f)  +  fjjCp,!.)  , 


(3-10) 


where 


tjCP,b  - 


(3-10a) 


and 


Zfi.v  -  -i 


for  £  e  I  (p) 


— —  2 

Ltt  zexp(-£  )  for  £  e  II  (p) 


,(D 


(3-10b) 


The  expression  (3-9)  for  T  may  then  be  divided  into  the 
four  terms: 


(1)  =  (1)  (1)  (1)  _  (1) 

1  1 II’I*  +  J I'll’  +  1 II’II’  1  I  II 


(3-11) 


with 


ra> 

IIT 


=  [dr  12  Jd612K12  I x-p 

CT  -*■  * 

erp<o 


(3-lla) 
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r^jj  ,  =  -4ttR‘ 


lVoofdr  12  Jd^1^2  fd612K12  ^1-2 
°  ti*M 


fj  (p-|lT12»  f  ii  (P“^iTi2*  ^2^  1  ' 


r(1) 

1  ini’ 


-4ttR  mn  go 

CO  Co, 


,(1) 

X  II 


DO 

JaT12  f^A2  fa&12K12  V 
0  ■+  *  _ 
ea-p<o 

T  £T)f^(p-£  T  £')  + 

Lrll  ^1T12,^1;III  Kp  n  12’^2; 

”  f  II  ^P”^iTi2,^P  f  II  ^p~^>  It12)  ^2  ^  ]  * 

« 

oo 

-47TR  mn^cD  /dr12  J  ^J2IC12  ^1*P 

?j-P<0 

[  f  I(0)  1T  H  •  *  1>  *  I?  <^1T  U  •  *  J>  * 


Note  that  the  condition  £^-p<;0  implies  that  in  (3-9] 


(3-llb) 


(3-llc) 


(3-lld) 


(3-12) 
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for  a  droplet  with  a  surface  whose  radius  of  curvature  is 
positive  everywhere.  This  has  been  used  in  obtaining 
(3-lld) . 

The  condition  ^*p<0  implies  that  molecule  1  is  imping¬ 
ing  upon  the  droplet  with  velocity  From  (3-8)  it  is 

evident  that  the  integral  vanishes  unless  molecule  1  suf¬ 
fered  a  collision  with  molecule  2  at  a  time  r 12  earlier. 

This  collision  is  characterized  by  the  perihelion  vector 
°12*  The  vel°cit^es  42  are  the  velocities  after  this 
collision  and  ^ vel°c^'ties  prior  to  this 
collision.  The  contribution  to  the  integrand  is  dependent 
on  whether,  prior  to  this  collision,  molecule  1  and/or 
molecule  2  originated  from  the  droplet  or  from  infinity. 

We  thus  consider  the  motion  of  two  isolated  molecules 
in  the  presence  of  the  droplet.  All  of  the  possible 
collision  sequences  involving  three  or  more  successive 
collisions  between  two  molecules  and  the  droplet  are 
schematically  represented  by  the  four  diagrams  of  Fig.  2. 

The  lines  indicate  the  trajectories  of  the  particles; 
the  circle  represents  the  droplet.  The  position  of  the 
molecules  at  time  r  =  0  is  indicated  by  the  dots;  the 

■f 

velocities  of  molecules  1  and  2  at  zero  time  are  and 
£2.  At  r  =  0  molecule  1  is  impinging  upon  the  droplet, 
and  in  all  diagrams  the  time  is  taken  to  increase  when 
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Pig.  2.  The  collision  sequences  among  two  vapor  molecules  and  the  droplet 
that  are  associated  with  the  coefficient  r  ^  of  the  first  inverse 
Knudsen  number  correction  term.  The  lines  represent  the  trajec¬ 
tories  of  molecules  1  and  2  and  the  circle  represents  the  droplet. 
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the  trajectories  are  traversed  in  the  direction  of  the 
arrows.  We  also  require,  in  all  of  the  collisional 
sequences,  that  molecules  1  and  2  collide  at  time  r  =  -t^2. 
The  diagrams  in  Pig.  2  differ  in  the  dynamical  history  of 
the  molecules  prior  to  the  collision  between  1  and  2.  In 
the  R-  collision  (recollision)  we  require  that  molecule  1 
originate  from  the  droplet,  in  the  C-  collision  (cyclic 
collision)  that  molecule  2  originate  from  the  droplet  and 
in  the  RC-  collision  •  (recollision-cyclic)  that  both  mole¬ 
cules  1  and  2  originate  from  the  droplet.  In  the  H- 
collision  (hypothetical  collision)  we  require  that  2 
originate  from  the  droplet,  but  here  we  identify  the  velo¬ 
cities  f  <t2 '  prior  to  the  (12)  collision  with  the  veloci- 
■+ 

ties  £2  after  the  (12)  collision.  We  shall  refer  to 

such  a  (12)  collision  as  a  non-interacting  collision.  In 
each  of  these  collision  sequences,  we  wish  to  remain  non¬ 
committal  about  the  possible  occurrence  of  collisions  that 
are  not  indicated  explicitly.  This  interpretation  of  the 
diagrams  is  for  convenience  and  has  no  effect  on  the 
exactness  of  our  calculation.  For  instance,  in  the  R- 
collision  molecule  2  may  or  may  not  have  originated  from 
the  droplet  prior  to  the  (12)  collision  and  molecule  2 
may  or  may  not  impinge  upon  the  droplet  after  the  (12) 
collision. 
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Each  diagram  in  Fig.  2  determines  a  region  of  inte- 
gration  for  the  variables  4r  £2,  r12,  0l2  such  that  the 
collision  sequence  indicated  takes  place.  From  (3-lla) 
we  see  that  the  integrand  of  vanishes  unless 

particle  1  originated  from  the  droplet  and  particle  2  from 
infinity.  If  we  now  replace  the  distribution  functions  by 
the  appropriate  Maxwellians  and  integrate  over  those 
regions  of  the  variables  4^  42,  t  ,  C12  associated  with 
the  recollision  sequence  of  Fig.  2,  we  make  an  error  since, 
in  so  doing,  we  have  also  implicitly  included  integration 
regions  associated  with  the  RC-  collision.  This  contribu¬ 
tion  therefore  needs  to  be  subtracted.  We  thus  obtain: 


(3-13) 

The  notation  f  indicates  that  the  integration  region  is 

(R> 
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limited  so  as  to  satisfy  the  conditions  for  a  recollision 
to  occur.  The  other  collision  integrals  in  (3-11)  can  be 
related  to  the  diagrams  of  Fig.  2  in  a  similar  way  to  yield: 


(1)  —22  /"*'*■  -*■  .  +  A 

rr  ir  -  -47T  r  mnAJd^ idTi2d^2d8i2x:i2 


J 

( C ) 


V exp  ■ 

/  2  iiZ  \ 

-  (?i2  +  2*  ) 

*  -exp  ■< 

-(q2+ 

* 

.T  2 

1  J 

L 

- 

-2  2 

+  47T  R  mn  m 

CO  00 


lt 


4  exp  ]-  (q2  +  ^-) 
*7 


idr12d^2d^12'c12  ^l’^ 

(RC) 

I’2 


-exp  i 


-cei2n2'2) 


(3-14) 


rn),ir  "  -^"2h2»v. 


fal dT12d^2d812*:12  ^l'P 

(RC) 


+2 

n 

~ r  e^P  H 

L  f!i!  +  ii!)1 

■  -exp  - 

-<q2+  e2,2)V 

[  v  T  T  yJ 

l  jJ 

(3-15) 


'5^ jt  =  -47T_2R2mn  CD 
I  ll  »  » 


fa  ldT12d^2d^12'c12  ^l’P 

(H) 

4  exp  j  (,i2 + ii!y 


lt 


-exp  1 


-(?i2+  e22) 


(3-16) 
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Upon  rearranging  the  terms  in  (3-13)  through  (3-16)  we 
conclude  that: 


r(1)  =  PR1)+  FC1) 


+  r 


(i) 

H 


+  r 


(l) 

RC 


(3-17) 


where  F^,  and  are  given  in  Table  I.  The 

collision  integrals  F^,  are  associated,  respec¬ 

tively,  with  the  R-  C-  and  H-  collision  sequences  which 
involve  three  successive  collisions  among  two  molecules  and 
the  droplet;  is  related  to  the  RC-  collision  sequence 

with  four  successive  collisions.  It  is  shown  in  Appendix  A 
that  these  collision  integrals  are  convergent,  so  that  the 
coefficient  of  the  term  linear  in  a  is  finite. 
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TABLE  1 


The  colllelen  integrals  far  the  coefficient  r 


(O 


af  the  contribution  proportional  to  a. 
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IV.  Second  Correction  Term  to  the  Free  Molecular 
Approximation 

We  consider  the  solutions  f^(p,£)  and  f^(p,|)  as 
the  first  two  terms  of  an  expansion  for  the  distribution 
funct ion: 

tCp.b  -  f (0)(p,b  +  f(1)(p,?>  +  f(2)(p,i>  +  ...  (4-D 


If  we  substitute  this  expansion  into  the  Boltzmann  equation 
(2-7)  and  neglect  terms  which  are  obviously  of  higher 
order,  we  obtain  the  equation: 

\  .  M<2)  Cp.1i)  -  a  fit  d6  k 
*1  ^  1  a  Ja^2a0l2IC12 


{[fc°) 

(1) 

(p,  1$ 

4-  f(1> 

(0) 

(P*?^  + 

-  f(0) 

Cp,  fx)f 

(1) 

(p»l2^ 

_  f(l) 

(p,ipf 

(0) 

Cp  *  ^  2^  J  + 

+  [f(0> 

Cp,t{)i 

(2) 

cp42) 

+  f(2) 

(2) 

(p,ip  + 

-  f(0) 

Cp.ljt 

(2) 

(p>f  2) 

_  fC2) 

(p.f 

(0) 

+ 

4. 

(P'tp* 

(1) 

Cp,t$ 

_  f(l) 

(Pjlpf 

(1) 

In  the  previous  section  we  showed  that  leads  to 

a  correction  term  linear  in  the  expansion  parameter  a.  If 
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f  2)  O 

we  now  were  to  assume  that  f  is  proportional  to  a  ,  we 

would  retain  only  the  terms  involving  on  the  right 

hand  side  of  (4-2) .  This  equation,  thus  truncated,  would 

yield  in  analogy  to  (3-6)  : 

t 

*W>  <?■«!>  ‘  a/dT12/^2d812K12 
0 

•  i*(0)(P-l1T12.?i,)f<l>(P-?lTl2-?2)  + 

-  ft0)^ri2-?i^(1)^-Vl2.?2>  + 

-  f<1)fi5-l1T12,?1)iCO)(?-I1T12,t2)]  ,  (4-3) 


where  the  upper  limit  t  of  the  t12-  integration  is  again 
given  by  (3-4) .  The  resulting  contribution  to  the  flux  r 
would  then  read: 


(0) 


+  ar 


(1) 


a2r<2> 


(4-4) 


with 
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r(2)  -  -Wm„„a.Ma-iy’d?1I1.p/dr12  /d^do^ c12 

ex-p<o  o 

■  [f(0)(H1T12,I1')l(1)(p-V12,|^)  + 

-  *<0)(p-I1T12,11)f(1)(p-11T12,i2)  + 

-  fw><p41T12,e1)f<0)(p-11T12,t2)J  .  u-5) 

In  the  previous  section  we  showed  that  is  deter¬ 

mined  by  integrals  associated  with  three  or  more  successive 
collisions  among  two  molecules  and  the  droplet.  The 
expression  (4-5)  can  similarly  be  analyzed  in  a  manner 
that  will  be  discussed  later  in  this  section.  It  turns 
out  that  (4-5)  can  be  decomposed  into  a  set  of  integrals 
associated  with  four  or  more  successive  collisions  among 
three  molecules  and  the  droplet.  The  leading  terms  are 
determined  by  the  integrals  corresponding  to  the  sequences 
of  four  collisions  depicted  in  Fig.  3.  However,  as  shown 
in  Appendix  A,  the  phase  space  associated  with  these  events 
behaves  as  for  large  values  of  r12  and  thus  the  inte¬ 
grals  diverge  logarithmically.  Our  original  assumption, 

'2) 

therefore,  that  1  “  yields  a  contribution  proportional  to 
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Fig.  3.  The  collision  sequences  among  three  vapor  mole¬ 
cules  and  the  droplet.  They  are  associated  with 

r(2:>  in  (4-5)  and  with  the  coefficient  of 

2 

the  correction  term  proportional  to  a  lna. 
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2 

a  ,  is  not  justified  and  we  need  to  consider  the  full 
equation  (4-2) . 

Since,  on  the  right  hand  side  of  (4-2),  the  Boltzmann 
collision  operator  operates  on  the  unknown  function 
solving  this  equation  poses  a  formidable  task.  It  is, 
however,  possible  to  extract  the  contribution  to  the  mass 
flux  which  is  of  lowest  order  in  <x.  For  this  purpose,  it 
suffices  to  retain  on  the  right  side  of  (4-2)  the  expres¬ 
sions  involving  and  the  term  which  is  directly 

proportional  to  f ^  (p,£j)  : 

a  jrd12d812x:12f(2)Cp,l1)f(0)(^)|2)  .  (4-6) 


We  thus  consider  the  equation: 


bp 


=  a  /<U2dS12«:l2  [f(0)(?,I1')f(1)(?,I2')+  f(1)  Cp,%P  + 


(4-7) 


where  v(^)  is  the  collision  frequency: 
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fd^2d^12IC12io>^2?  *  (4-8) 

In  writing  (4-8)  we  have  approximated  f^°^(p,|2>  in  (4-6)  by 
the  Maxwell  distribution  fM(£2)  •  The  justification  for  this 
will  be  explained  shortly. 

The  simplified  equation  (4-7)  can  be  readily  integrated 
to  yield: 

t 

(o')  ->  ■+■  r  r  a  —  q,vti2 

f  (P’ii)  “  a  J  dTi2  J  d^2dc12*r12e 
0 

+  f ( 1  *  ^ -i  ^  12 ,  i  p  1 (o)  <p-i  lT  12 ,1 2')  + 

-  f(1>(p-i1T12,l1)f(0)  (P-11t12,I2)]  .  (4-9) 

This  result  is  the  same  as  (4-3)  except  for  the  inclusion 
of  the  exponential  factor  e"avTl2.  This  damping  factor 
will  insure  that  the  resulting  collision  integrals  are 
convergent . 

With  this  damping  factor  introduced  into  (4-5) ,  we 
conclude  that  for  large  values  of  r ^  the  leading  term 
in  the  second  correction  to  the  mass  flux  r  will  vary 
as: 
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a*f' 


-avr 


dT12e 


12  _A 

T 


12 


Aa  lna  +  terms  of  higher  order  in  a. 

(4-10) 


The  situation  thus  appears  to  be  quite  analogous  to  the 

exponential  damping  of  the  collision  integrals  in  the  theory 

of  the  density  dependence  for  the  transport  coefficients  of 

a  gas  [19,20]  .  While  the  expression  (4-9)  for  f^  does  not 

represent  a  complete  solution  of  the  original  equation  (4-2), 

it  does  contain  all  necessary  information  to  determine  the 

coefficient  of  the  contribution  to  the  mass  flux  proportional 
o 

to  a  lna.  The  reason  for  this  is  twolold.  First,  the  terms 
containing  give  a  contribution  to  r  which  is  of 

3 

order  a  .  Secondly,  the  neglected  terms  in  (4-2)  involving 
f^  and  also  the  replacement  of  f^(p,f2)  in  (4-6)  by 
f„(e2)  modify  the  details  of  the  exponential  damping  factor, 
but  do  not  change  the  coefficient  A  in  (4-10) .  The  mathe¬ 
matical  analysis  to  substantiate  these  assertions  is  quite 
complex  and  is,  to  the  extent  feasible,  presented  in 
Appendix  B.  In  the  main  body  of  this  report  we  shall  con¬ 
centrate  oil  formulating  the  collision  integrals  that  will 

enable  us  to  calculate  the  coefficient  of  the  contribution 

2 

proportional  to  a  lna. 

We  thus  consider: 


(0)  +  r(l)a  +  r(2)a2lna  + 


where : 


(4-11) 
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f(2)-4IR2mnAli.  /d?  Ji'P  /dT12/d«2d512':12e  .  ”  “ 

t1-P< 0  ° 

♦  f(1)(p-i1T12,i1')J(0)(p41T12,1^)  + 

-  *(0)(?-i1T12.I1)*(1)<p-e1T12,iil)  + 

-  ] ,  (4-12) 


is  obtained  by  substituting  (4-9)  into  (2-10) .  This  result 
differs  from  (4-5)  due  to  the  inclusion  of  the  exponential 
damping  factor. 

In  order  to  avoid  the  repetition  of  unnecessary 
details,  we  introduce  a  shorthand  notation: 

00 

[ f  (0>  (1) f  (1)  (j)]  - -Wmn^  lim  ^ /<!?& ■  pf Jt12  /d|2d S12 

li-frc*  ° 

•  ici2eavTl2f<0)(^ri2'5i)f<1)(p-^ri2’^j)  1  (4"13) 

so  that 


f<2)- 

f  (°)(r)f  (D  (2.)  j 

+ 

[f  (1)(l?)f  (0)  (2  r)  J 

- 

f  (l)f  (1)  (2)  j 

- 

f  (!)  (1)f  (°)  (2)  J 
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where  all  distribution  functions  in  the  integrands  are  to 

A  -f 

be  evaluated  at  the  position  p  -  p-^T^. 

~  (21 

We  shall  demonstrate  that  this  coefficient  r  is 
determined  by  a  set  of  integrals  associated  with  sequences 
of  collisions  among  three  molecules  and  the  droplet.  For 


this  purpose  we  first  write,  in  accordance  with  (2-14)  and 


(3-6a)  f 


and  f(1>  =  f<»  + 


to  obtain: 


r(2)  -  [fj0)  (l’)f jX)  (2’)]  +  [fj15  (l*)fj0)  (2’)]  + 

+  [f  J0)(l»)f  JJ)(2t)]  +  [f{1)  (l’)f  (21)]  + 

+  [f£^  (l')**1*  (2*)]  +  [f^  (l*)f<0)  (21)]  + 

+  [fJJ)(l')f^)(2’)]  +  [f<l)(l,)f<°)(2,>]  + 

-  [f^0)(l)  fj1)(2)  ]  +  f  j0)  (2)  J  + 

-  [fj0)(l)  f^(2)  ]  -  [fj15(l)  fjj)(2) 

(4-15) 

Here  we  have  used  the  fact  that  in  the  integrand  of  (4-12) 

?!  e  Klp-I^jgl)  so  that  f (0)  (P-?1T12»?1)  “  f  j0)  (p-?iTi2»1i) 
and  f  (1)  (p-?1T12,?1)  -  f  j15  Cp-%iri2^i>  ■ 

As  a  second  step  we  decompose  f^1^,  given  by  (3-6b)  in 
the  same  manner  as  was  done  in  (3-11)  for  .  We  then 

obtain: 
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+  (i)  -  f 

+  1IJII,II,V  '  II,I  IIU;  ’  (4-16) 


-r 

where  we  have,  lor  £ ^  e  I  ( | p-^r^ I ) " 


aJdTi3  Jd^3dSi3'ci3 


(i) 

tI, II 1 1 ’ UJ 


0 


■  rf<°>  rp_i  T  T  ■?  u\ f  (°)  T  r  i ")  + 

L  II  n  12  M  i3’M  }TI  Kp  n  12  13'^3  ;  + 

-  fIl(^lT12-1lTi3-«l")*I(^lT12-?lTi3^3">]  *  (4-16a) 

00 

fI,I,II,(l)  "  a  fdT±3  Jd^3d^i3K±3 


fO\  *  •*  -*■  -*■  (n'i  „  ->■  -*■  ■*• 

,fI  (p“^lT12~^iTi3’^i^fII  (p_^lT12“^lTi3'^3 + 


fI(p_^lT12"^iTi3’^i  ")fII(p"^lT12"’^irri3,^3")  ]  '  (4  16b) 


u) 

fI,II,II’^l)  “  a  fdTi3  Jd^3d^i3Ki3 

0 

[f  jl^  (P“llTl2“^iTi3»^iM)f (p-?iTi2"^iTi3,^3,,)  + 
flI(p-^1Tl2-^iTi3^i")fiifp-€1T12_^iTi3>^3M)  ]  »  (4-16c) 
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TO 

fl”l  II(i)  =  a/dT13 
0 

[fI<0>  <^lT12-^iTi3>4i»l(?)  (P-Vl2-Vi3’V  + 
fI^p"^lT12~^iTi3’^i'fII^p_^lT12_^iTi3’^3^  ]  '  (4-16d) 


"T 

For  £.  e  II  (|  p-£jT12|  )  we  take,  just  as  in  (3-6b) : 


ftt>  (i) 

*if  in* U) 


rl,  I'  II 1  u; 


f  (1)  (jn 

1,11'ir 


Qj  a  o 
*1,1  IIK  }  u- 


(4-16e) 


In  (4-16a)  through  (4-16d) r  o..^  is  the  perihelion  vector  of 
a  collision  between  molecules  i  and  3  taking  place  at  the 
position  p  =  P“?iTi2‘”^iTi3 '  The  initial  velocities  prior 
to  this  collision  are  and  £3"  and  the  final  velocities 
after  this  collision  are  £ .  and 

1  O 

Similarly  the  function  f  »  given  by 

(3-6c) ,  is  decomposed  as: 


f 


(1) 

II 


(p-ilTl2>£i) 


;(1) 

ii,  in 


,  (i)  +  f 


(i) 

ii, i'ii 


(i)  + 


+  f 


(1) 

II, II'II 


(i) 


f (i)  _ 

II, II  Iw 


,(1) 

II, II  II 


(i) 


(4-17) 
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where  for  ^  e  1 1( | P~41T12 P 

T(P“t 


(i) 

*11, II'I* 1  ; 


“/  dTlzJdhdZi3Ki3 


fII(p-11T12“^iTi3'^i")fiCp-1iTi2  ^iTi3,^3,) 


(i) 

1 1 1 , 1 '  II’  u; 


T(p  ^1T12* 


aJ  dr  i3  /d^gdS  i3^i3 

0 


*  [fI  ^  ^  ^lT12-^iTi3,^i")  flV  (^"^lrl2  ^iTi3’^3  ^ 

fl(p-ilT12“IiTi3^i,')fii(p“IiTi2"^iTi3»^3")  ]  ’ 


f  (i) 


T<p-iiT12,ii> 

“/  dTi3jdhd613Ki 


[«ff  <P^lT12^iTi3'?i"),I<?>  <e-Vl2-Vi3'V> 


-  fu<P"?lTX2"^lT13’^i  "':fII^"^lT12-^iTi3’^3’') 


(4-17a) 


(4-17b) 


(4- 17c) 
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(1) 


T<p-?r12,ii) 

(!)  =  aj  dTi3  y'd^3d°i3lci3 


•f  A 


*11,11  I 

[fII^  ^“^lT12“^iTi3’^i^fI  ^  ^“^lT12'^iTi3’^3^  + 
f  II(^lT12~^iTi3J  ^  fI^"^iT12“^iTi3J^3)  ]  *  (4 


,(1) 

II, II  II 


i3  J^3'4Vi3"’i3 


(i)  -  a/  dri3  y j^0ddJ0K 


[f  (p-€1T12*^iTi3*^i)fii^  (p-?1T12-IiTi3»I3)  + 
fII^-^lTl2“^iTi3,^i)fII^“^lT12"^iri3,^3>  ]  * 


Just  as  in  (3-6c) ,  it  is  again  understood  that  for 

?.  c  KIp-I^P 


fn!irr(1»  "  fii! j’ii’  (i)  “  fii! ii'ii’ (i> 


f  (i)  =  f  =  0 

II, II  Iv  ;  ZII,II  IIKX) 


In  order  to  investigate  the  structure  of  the  colli 

~  (2) 

sion  integrals  determining  r  ,  we  consider  as  an 
example  the  term: 


-17d) 


-17e) 
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f(°)  (2') 


2  1 
-4ttR  mn  a)  lim  5— — 
m  «  _  ma 

a-*-0 


*  ■fihh-P  ,/dT12  Mid&12K12e 

e1-p<«  0 


-avT12  ? 


r. 


|^"oi  q  /  t  q^o  t 


d' 


2  *  2  j  j  2  1 3’'*  2  f  ^ 


fI  ^(P“£iTi2>?p  [  fII  ^  ^1T12“^2T2,3’^2  1T12_^2T2,3 ’^3 

tII^~^lT12~^2T2'3’^2 fI^-^lr12~^2T2'3^3  ]  ’ 


(4-18) 

which  is  obtained  when  (4-16a)  is  substituted  into  the  first 
terra  of  (4-15) .  The  integrand  in  (4-18)  corresponds  to 
those  phases  of  molecule  1,  2  and  3  such  that: 

(a)  at  p  -  p,  corresponding  to  the  time  r  »  0, 
molecule  1  is  impinging  upon  the  droplet; 

(b)  at  p  -  p— 4-jT corresponding  to  the  time  T  =  ~T \2* 
molecule  1  collides  with  molecule  2;  (the  velo¬ 
cities  of  these  molecules  prior  to  this  collision 
are  and  £2’  and  after  this  collision  ^  and  |2)  ; 

(c)  at  p  »  P~%lrl2~^2  t2'3’  corresponding  to  the  time 
t  -  “(Ti2  +  t2'3^  molecule  2  collides  with  mole¬ 
cule  3.  The  velocities  of  these  molecules  prior 


to  the  collision  are  and  !*„"  and  after  the 

2 '  and^3 


■2  - b3 

collision  are  |9’  and  ; 
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(d)  before  the  collision  between  molecules  2  and  3, 
the  velocity  of  2  must  be  within  region  II,  which 
means  that  molecule  2  originated  from  the  droplet. 
This  collision  sequence  is  shown  schematically  in  Fig.  4 
using  the  diagrammatic  notation  developed  in  the  previous 
section.  The  integral  (4-18)  represents  an  integral  over 
all  phase  points  of  molecules  1,  2  and  3  which  allow  for 
the  occurrence  of  the  sequence  depicted  in  this  diagram. 

In  the  previous  section  we  referred  to  the  collision 
sequences  that  contributed  to  the  first  coefficient 
as  R- ,  C-,  H-  and  RC-  collisions.  In  order  to  specify  the 
collision  integrals  for  the  coefficient  r  we  need  a  more 
systematic  notation.  We  shall  indicate  an  interacting 
collision  between  molecules  i  and  j  by  (ij)  and  a  non¬ 
interacting  collision  between  i  and  j  by  (ij)  .  A  collision 
between  molecule  k  and  the  droplet  will  be  denoted  by  (Ok) 
where  O  refers  to  the  droplet.  A  collision  sequence  is 
represent  by  a  lef t-to-r ight  juxtaposition  of  these  sym¬ 
bols.  In  this  notation  the  R-,  C-  and  H-  collisions 
introduced  in  the  previous  section  can  be  represented  by 
(01)  (12)  (01),  (02)  (12)  (01)  and  (02)  (12)  (01) ,  respectively, 
and  the  diagram  of  Fig.  4  by  (02) (23) (12)  (01) .  Thus  the 
contribution  (4-18)  to  r  becomes: 
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Fig.  4.  Example  of  a  sequence  of  four  collisions  among  three  vapor  molecules 
and  the  droplet.  In  this  example  molecule  2,  which  has  originated 
from  the  droplet,  first  collides  with  3  and  subsequently  with  1  in 
such  a  manner  that  molecule  1  will  impinge  upon  the  droplet. 
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f(2) 

1  (02)  (23)  (12)  (01) 


f<0)  (!•)*<»,, I,  (2')^ 


=  ~47rR2mnA  lira  ^  /df idT12dI2d3l2dT2, gd^dOg , 

/rx/N  \  y  y.  /A.x 


(02)  (23)  (12)  (01) 


-v  ^  “ClVT  22  3  „  r  o3  ^ 

K12'c2,3^1*^e  7r'Texp{-e]'2}[n*7r"3T*_Texp{-(-|-+  £3"2)} 


-  7r-3exp{-  (^  ,,2+  £  ”2)} 


(4-19) 


where  the  symbol  j  means  that  the  integration 

(02)  (23)  (12)  (01) 

extends  over  all  phases  of  molecules  1,  2  and  3  corre¬ 
sponding  to  the  collision  sequence  (02) (23) (12) (01) . 

If  we  substitute  (4-16)  and  (4-17)  into  (4-15)  we  find 
~  (2) 

that  r  appears  to  be  a  sum  of  fifty-three  terms  each  of 
which  may  be  associated  with  a  sequence  of  collisions 
among  three  vapor  molecules  and  the  droplet.  However, 
many  of  these  terms  correspond  to  sequences  that  involve 
more  than  four  successive  collisions.  As  an  illustration 
we  examine  the  two  terms  [  f  jj  '  d  ’)f  (2')]  and 

[f j jj i j i (1 ' ) f (2 ’ ) ]  which  correspond  to  the  collision 
sequences  shown  in  Fig.  5.  These  two  terras  both  require 
a  sequence  of  collisions  which  in  our  notation  can  be 
indicated  by  (01)  (13)  (12)  (01)  .  The  term  [fj  n  »  jj.d')  f  j°*(2  ')] 
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Fig.  5.  Examples  of  sequences  of  five  collisions  among  three  vapor  molecules 
and  the  droplet. 

(a)  Diagram  corresponding  to  [  f  j1  jj.  ,  ,  (1 ' )  f  (2' )  3  , 

(b)  Diagram  corresponding  to  [f  (l')f  ^  (2')1  . 
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requires,  in  addition,  that  £3"  e  II  Cl  P“I1T12~^ i  ,Ti 1  3 I  )  which 
means  that  not  only  1,  but  also  3  must  have  originated  ini¬ 
tially  from  the  droplet.  The  term  [f  ^  ^  jj  ,  j. ,  (1  * )  f  (2‘  )  ] 
does  not  involve  the  latter  condition,  but  requires  instead 
e  II  (Ip-^T^I)  this  constraint  implies  that  the  phase 
of  molecule  1  must  be  such  that  it  also  would  have  origi¬ 
nated  from  the  droplet  if  we  ignore  the  velocity  change 
caused  by  the  (13)  collision.  Both  terms  therefore  involve 

constraints  related  to  five  collisions.  In  fact,  it  can 

~(o\ 

be  readily  verified  that  terms  in  r  which  contain  either 

fj®^  and/or  any  of  the  following: 

*(1)  -(1)  -(1) 

^  II,  I '  II '  ’  fII,II,II”  fII,II  II’  a11  corresP°nd  to 
sequences  of  five  or  more  collisions  among  three  molecules 

and  the  droplet.  However  as  discussed  in  Appendix  A  such 

terms  lead  to  integrals  of  the  form: 


lim 

cr*-0 


TSa  J't 


-avr 


dT12G 


12 


12 


2+y 


(4-20) 


with  y  2-  0  and,  therefore,  do  not  yield  a  contribution  to 

~  (2)  2 

r  In  order  to  determine  the  coefficient  of  the  a  lna 

term  it  is  sufficient  to  retain  only  the  collision  integrals 

related  to  four  successive  collisions.  We  thus  obtain: 
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r  (2)  =p  (2)  +f(2>  +p(2)  _  . 

r  1  (02)  (23)  (12)  (01)  +  1  (03)  (23)  (12)  (01)  +  1  (03)  (23)  (12)  (01)  + 

■  P  (2)  ,p(2)  ,  p  (2)  + 

+r(01)  (13)  (12)  (01)  +1  (03)  (13)  (12)  (01)  +  1  (03)  (1 2)  (12)  (01) 

(2)  .  « (2)  ,=(2)  _ 

+r(02)  (23)  (12)  (01)  +  1  (01)  (13)  (12)  (01)  +  1  (02)  (23)  (T2)  (01)  + 


■  P  (2)  _ 

+I  (03)  (23)  (12)  (01) 


?<2)  —  — 

1  (03)  (23)  (12)  (01) 


.  pf2)  _ 

+  1  (01)  (13)  (12)  (01) 


+ 


+r 


£  (2) 


(03)  (13)  (12)  (01) 


^  =  (2>  _  .=(2) 

+  "  (03)  (T2)  (12)  (01)  +  r  (02)  (22)  (T2)  (01) 


(4-21) 


where 


p(2) 

1  (02) (23)  (12) (01) 

p(2) 

1  (03)  (23)  (12)  (01) 

P  (2)  _ 

1  (03)  (23)  (12)  (01) 


P<2) 

1  (01) (13) (12)  (01) 

p(2) 

1  (03) (13)  (12)  (01) 

?(2) 

3  (03)  (13)  (12)  (01) 
?(2)  ^ 

r  (02)  (22)  (12)  (01) 
p(2) 

r(01)  (T2)  (12)  (01) 

1  (02) (23)  (12)  (01) 

f(2)  _ 

J  (03)  (23)  (12)  (01) 


-  +(f 

-  +[* 


-  “I* 


=  +lf 
=  +[f 


=  -If 


-  -If 


-  -[f 


-  -If 


"  ~[f 


9 

(4-21a) 

9 

(4-2  lb) 

ro')'!:1,  ii»'m 

9 

(4-21c) 

f 

(4-2 Id) 

i^rir  (1,)fi(0)  (2,)1 

9 

(4-21e) 

1*1  jjd*)^05  (2 1 ) ] 

9 

(4~21f) 

i0><i,’iu!ii  i(z'» 

9 

(4-21g) 

fiiii 

i 

(4-2 lh) 

I0) 

t 

(4-211) 

* 

(4-21  j) 
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p  (2) 

1  (03)  (23)  (12)  (01) 

-  +[f^0)(l)f^)I  IX(2)] 

(4-21k) 

-v  {2) 

r(01)(13)(12)  (01) 

=  ■Ifi!iri,(1):fiC0)<2)1  * 

(4-21£) 

P<2) 

*  (03)  (13)  (12)  (01) 

_[fi, r ii* (1)f i°* ’ 

(4-2  lm) 

~  (2) 

r  (03)  (13)  (12)  (01) 

-  +[f^)I  u(l)f{0)  (2)) 

(4- 2  In) 

p  (2) 

1  (02)  (23)  (12)  (01) 

-  +[f1(0)  (l)f^|n  x  (2)  ]  , 

(4-21o) 

~  (a} 

We  conclude  that  the  coefficient  r  is  determined  by 
a  sum  of  fifteen  collision  integrals.  They  are  associated 
with  the  diagrams  shown  in  Fig.  3  which  represent  all  of 
the  possible  sequences  of  four  successive  collisions  among 
three  vapor  molecules  and  the  droplet.  The  explicit 
expressions  for  these  collision  integrals  are  presented 
in  Table  II.  We  note  that  the  signs  of  the  collision 
integrals  depend  on  whether  the  number  of  non-interacting 
collisions  is  even  or  odd. 
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TABIC  XI  The  collision  integral* 


e  :  +4«”  ^  R2!!!^ 


r(omzj)(12XOl>  ■  -a  u"  lb  /‘-aVuVs  ^-.‘S  “M-C^l2  ♦  «;  » 

°*  («2)(23>(12)(01) 

8  (03) (23) (12) (01)  ■  *6  lb  /‘‘Va'ia’Va  V"  +  <1  >1  tj’2) l] 

‘r>0  (03)  (23)  (12)<«1) 

{2  a 

r(03)  (23)  (12) (01)  *  +S  lla  lb  /<h,2,3,la';,3  V*  •*pt-Ca«Tlz  +  t|  ))  [55777  «*p(-(«J!  +  5I  >>  -  «jM(i  +  t]»J 

a"°  (03)  03)  (12)  (01) 

2  r  2 

r(OX)(l3)(12)C01)  •  "fi  ^  lh  /d*  1,3<12,CI,3  VS  +  C^)  I  [y5t7T  +  c”’))  -  «»<-(("  *  c"’»] 

r(Q3)  (13)  (12)  (01)  *  -«  u"  lb  /ta  1V.2V1  «•  S  «»P<-<“VTJ2  +  t;S>)  [^l7I-«P<-(t;,1+1T  >1  -  «P<-(t"2+ 

(03) (13) (12) (01) 

(i2) (01)  -  +*  ih  /d* i'j*i>Vi  V5  “pM™,,  +  ^’>1  [55177  “pMeJ*  -  -xpt-te;2  +  Cj)>] 

“r’°  (03)(23)(12) (01) 

t>2 

(121(01)  ■  +a  lb  f  ** 2V1JV1  V®  «*p«*c«,2  +  t|2)>  [55777  «p:-(5?*  c*))  -  «pi-tc’a  +  ti»] 

a"  (02) (23) (12) (01) 

t’2 

(12)  (01)  ■  +a  U“  lb  /  dl*| 1 3d|3*l  *3  V*  **p{'<“*rl2  +  {i2))  [55177“?  -«t5‘+  C3)}  *  +  t!»] 

(01) (13) (12) (01) 

c”2  a  2 

U2) <oi)  *  +s  Um  lb  /^ai'iafai  V*  €*»(-6,v’r,j  +  05>  [tSt77  “P1"^  +  «i‘  »  "  “?<-<«!'  +  5 >] 

°~°  (02) (23) (12) (01) 

^m2 

(72)(01)  "  lno  f  Jlj V}  «xp{“(<*VT12  +  Jt*37T  ^  «p{-(C2  +  C,  ))J 

(03)  (23)  (12)  (01) 

(12)(0J)  *  '*  “!  ih  /dl,a3llJ^21  V®  ♦  C2))  [f?77r  “P<-«!  +Tf)>  -  «pl-tt|  +  tj)>] 

(03) (23) (12)(01) 

t 1  1  2 

(i2)  (01)  ■ +a  “**  ih  /  dxi3‘u^j  v®  ■«»(-<««„  +  k2>)  [55777  “p<"<tt  +  «rJ>}  ■  “p<-«rl  + 

3  (01) (13)(12) (01) 

(IT)  (01)  *  +s  “5  lib  /  d*:i‘ia^5  VS  «*»f-CovT13  +  Ea)l  [55777  *xpl-(tJ,2+5r  >)  -  «xp(-(t”J+  c”J))l 
(03) (13) (12) (01) 

(TT)  (01)  ’  -*  11"  /  d*i3':a,cii  «i-»  •»?«-<■«„  +  t2)>  [55772  «Pt‘t^  +  57  »  -  *»pt-U*  +  t2,)»] 

a*0  (03)  (H)  02)  (01) 


/  4«i3  i/i<idriadV5iad,i3',Vd13 


n2 

iJ  r  .»  ^a  ..a  • • a  1 . a  7 

+  >>  [55177  “Pl-fj*  +  t»  »  -  «pMtj  +  i,  »J 

2  l|2 

‘*v,ia  +  {!  11  [r5777  -“Pf-tt,  ‘  )>  -  «»pt-(tj,2+  2J  l] 


r (03) (13)  (12)  (01)  ’  44  lia  /d*l’ 3*1 1*1*3  *1*»  **P(-<a'lT1j  *  t 

(03)<13> (12) (01) 


r«2)(23>(12)(01)  ‘  +(1  “■  lb  17*2*1  *1*  “PMo^,  * 

(02) (23) (12) (01) 

r(01)(13)(12)(01)  ‘  +a  UJ|  lb  ■^d*ll  3d, 2*1*1  *“P<*<a»r,a  + 

(01) (13) (12) (01) 

r (02) (23) (T5) (01)  *  +s  u=  lb  ftx2}‘i7*_1,  «1-P  “P(-(av*,a  +  «! 

(02) (23) (12) (01) 

r(03)(23)(T2)(01)  ’  44  la!  lb  /  dlJ  1'  e*P(-tavTij  +  * 

(03X23)  (12)  (01) 

r(D3)(TI)(T2)(oj)  *  'a  “J  lb  f  d*a3^u]|2i  ^1'P  *”P{-<a"Iia  +  t 

(03) (23) (12)(01) 

r(01)(l3)(12)(01)  "  +S  ^  Tnl  /dxl3'ls^3  *1'P  a*l>(-<a»t1,  +  0 

(01)  (13X12)  (01) 

r(03)(13)(l2)(01)  *  +S  U°  lb  /dx:i‘l2^3  'rS  «»f-CavT13+  E 

(03)  (13X12)  (01) 


rE03)(l3)  (12)  (01)  "  'S  J1"  lno  /dxn':a‘:  11  V°  ax?1'(avT,j 
a*°  (03)  (H)  02)  (01) 

r{o2) (23) (72X01)  ’  -8  u"  lb  /dxJj:i?^  V»  •*P<-ta'-’n 

(02) (23) (12) (01) 


E1 

+  c!n  [t?T7T  *xPf"l55  +  -  «pt-(Cj  +  c“>)] 
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V.  Discussion 

In  this  report  we  have  considered  an  expansion  for  the 
mass  flux  r  in  the  nearly  free  molecular  regime  having 
the  form: 


(0) 


+  r(1)a 


+  r(2)a2lnct  + 


(5-1) 


where  a  is  the  inverse  Knudsen  number.  The  coefficient 

T(1)  is  determined  by  a  set  of  integrals  associated  with 

sequences  of  collisions  among  two  vapor  molecules  and  the 

~  (2 ^ 

droplet.  The  coefficient  f  '  of  the  logarithmic  term  is 
determined  by  integrals  associated  with  the  asymptotic 
behavior  of  collision  sequences  among  three  vapor  mole¬ 
cules  and  the  droplet. 

The  structure  of  the  expansion  (5-1)  suggests  an 
analogy  with  the  density  dependence  of  the  transport 
properties  of  a  moderately  dense  gas.  As  discussed  in  an 
earlier  technical  report  [8],  a  transport  coefficient  of 
a  gas,  such  as  the  viscosity  r),  should  be  written  as: 

2 

7)  =  T)0  +  T^n  +  7)2n  In  n  +  .  .  .  ,  (5-2) 

where  n  is  the  number  density.  The  value  r)0  of  the 
viscosity  in  the  low  density  limit  is  given  by  the 
Chapman-Enskog  theory  as  an  integral  over  the  parameters 
of  a  collision  between  two  molecules  [13] .  The  coefficient 
of  the  first  density  correction  involves  the  effect  of 
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correlations  in  the  position  and  velocity  variables  of 
three  molecules.  These  correlations  are  of  both  a 
statistical  and  a  dynamical  nature.  The  statistical 
correlations  refer  to  correlations  in  configuration  space 
irrespective  of  the  velocities  of  the  particles;  for  a 
gas  of  hard  spheres  they  reduce  to  excluded  volume  contri¬ 
butions.  The  dynamical  correlations  are  brought  about  by 
sequences  of  successive  collisions;  their  contribution  to 
is  given  by  collision  integrals  associated  with  three 
and  four  successive  collisions  among  three  molecules  [10] . 

The  sequences  of  three  successive  collisions  that 
appear  in  the  calculation  of  the  first  correction  term  rjj 
are  represented  schematically  in  Fig.  6.  In  earlier 
reports  [8,10] ,  we  have  referred  to  these  sequences  as  a 
recollision,  a  cyclic  collision  and  a  hypothetical  colli¬ 
sion.  The  events  indicated  by  R,  C  and  H  in  Fig.  2  are 
obtained  from  those  shown  in  Fig.  6  when  molecule  0  is 
replaced  by  the  droplet.  It  appears,  therefore,  that  a 
close  correspondence  exists  between  the  two  expansions 
(5-1)  and  (5-2) .  The  collision  sequences  to  be  consid¬ 
ered  in  the  evaluation  of  the  coefficients  of  (5-1)  are 
identical  to  those  previously  considered  in  the  density 
expansion  for  the  transport  coef f icients . if  we  identify 
one  of  the  molecules  with  the  droplet.  The  expansion 
parameter  in  (5-2)  is  actually  nO  =  C/(nC  )  which  is 
the  ratio  of  the  size  6  of  the  molecules  to  the  mean  free 
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Fig.  6.  Sequences  of  three  successive  collisions  among  three  molecules 
that  are  associated  with  the  first  density  corrections  to  the 
transport  properties  of  a  gas  [10] .  The  lines  represent 
schematically  the  trajectories  of  molecules  0,  1  and  2. 
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path  £;  the  expansion  parameter  in  (5-1)  is  the  ratio  of 
the  size  R  of  the  droplet  to  the  mean  free  path  i. 

This  correspondence  between  the  density  expansion  of 
the  transport  properties  of  a  gas  and  the  inverse  Knudsen 
number  expansion  of  the  properties  encountered  in  rarefied 
gas  dynamics  was  noted  earlier  by  Dorfraan  and  Sengers  [11] . 
The  density  expansion  (5-2)  for  the  transport  properties 
is  obtained  by  solving  the  generalized  Boltzmann  equation. 
Dorfman  etfll.  have  developed  from  this  same  equation  a 
theoretical  formulation  of  the  aerodynamic  forces  on  an 
object  in  a  rarefied  gas  stream  by  treating  the  object  as 
a  heavy  particle  [11,12].  Since  in  our  droplet  growth 
problem  the  number  of  molecules  in  the  vapor  phase  is  not 
conserved,  we  derived  the  mass  flux  r  by  solving  the  in¬ 
homogeneous  Boltzmann  equation  subject  to  the  appropriate 
boundary  conditions.  While  the  details  of  the  collision 
integrals  are  different,  it  turns  out  that  the  collision 
sequences  we  obtained  are  precisely  the  same  as  those 
derived  for  the  drag  coefficient  of  an  object. 

Collision  integrals  determining  the  first  inverse 
Knudsen  number  correction  to  the  drag  coefficient  of  a 
sphere  in  a  gas  stream  have  been  calculated  by  Kuperman 
and  Sengers  [21]  .  The  results  of  these  calculations  will 
be  reported  in  a  future  technical  report.  It  appears 
that  this  procedure,  with  minor  modifications,  can  be 
adapted  to  calculate  the  coefficient  in  (5-1). 
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APPENDIX  A 

PHASE  SPACE  ESTIMATES  FOR  COLLISION  INTEGRALS 

In  the  main  text  of  the  report  we  used  the  fact  that 
the  collision  integrals  determining  the  coefficient  are 

convergent  and  that  the  phase  space  volumes  associated  with 
the  collision  sequences  shown  in  Fig.  3  diverge  logarithmi¬ 
cally.  Here  we  discuss  the  evidence  for  these  assertions. 

The  situation  appears  to  be  completely  analogous  to  the 
phase  space  volumes  of  the  collision  integrals  for  the 
transport  properties  of  gases,  which  were  analyzed  by 
Dorfman  and  Cohen  [16] .  In  order  not  to  prejudice  the 
result  by  this  analogy,  however,  we  have  made  an  independent 
study  of  the  collision  integrals  derived  in  this  report.  A 
similar  analysis  of  the  collision  integrals  for  the  Knudsen 
number  dependence  of  the  force  on  a  macroscopic  object  in  a 
rarefied  gas  was  made  by  McClure  and  Dorfman  [12] . 

In  order  to  make  an  estimate  of  the  phase  space  vol¬ 
umes  associated  with  the  collision  integrals  for  we 

consider  as  an  example  the  recollision  contribution: 

00 

rRi;  -  fiT12 

o  (R) 

*  £  '2 

•  S-r  exP{-<-!-+  e2'2)}-  expl-di,^  e2'2)}]  .  (A-l) 

T*r  t 

■+ 

Since  the  integrals  over  the  velocities  ^  and  ^  are  well 
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behaved  as  a  result  of  the  presence  of  the  exponential 
factors,  the  convergence  of  depends  upon  the  behavior 

of  the  and  6^“ integrations  for  large  values  of 

It  is  convenient  to  change  the  sign  of  the  velocities 

and  i2  ar*d  t0  consider  the  recollision  event  of  Fig.  2 
with  the  direction  of  motion  of  molecules  1.  and  2  reversed. 
Let  us  choose  a  coordinate  system  in  which  molecule  1  is 
at  rest  at  the  origin  just  prior  to  the  (12)  collision  and 
in  which  the  Z-axis  is  taken  in  the  direction  of  -|^  (see 
Fig.  7) .  The  droplet  is  then  moving  in  the  Z-direction  and, 
at  the  time  of  the  (12)  collision  (t=t12) >  its  center  is 

located  at  r  -  ~P~^>iTi2  "  -^1T12  for  lar£e  values  of  t12. 
After  the  (12)  collision,  molecule  1  moves  in  this  coordi¬ 
nate  system  in  the  direction  of  with  a  velocity 
l-pi  s  ^l”^l-  In  order  for  *  to  catch  up  with  the  droplet, 
the  perihelion  vector  O^g  must  be  restricted  to  a  solid 
angle  proportional  to  r  -  |4jT12[  •  Thus  the  phase 

space  volume  behaves  for  large  values  of  t12  as: 


fdg 

12  J  QC 

(R) 


12 


(A-2) 


which  is  clearly  convergent  when  the  upper  limit  of  integra¬ 
tion  goes  to  infinity. 

The  integrals  and  also  involve  three  succes¬ 

sive  collisions  and  exhibit  the  same  asymptotic  behavior 
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as  r^15 .  The  collision  integral  requires  the  occur¬ 

rence  of  an  additional-  collision  and,  therefore,  converges 
even  more  rapidly.  In  fact,  it  turns  out  that  the  RC- 
collision  sequence  is  only  possible  when  the  (12)  collision 
takes  place  close  to  the  droplet. 

We  next  evaluate  the  phase  space  volumes  associated 
with  the  sequences  of  four  collisions  among  three  vapor 
molecules  and  the  droplet  shown  in  Fig.  3.  As  an  example 
we  consider  the  term: 


r(2) 

1  (02)  (23)  (12)  (01) 


00 

dT2'3  /dM^2d^3dS12 
(02)  (23)  (12)  (01) 


:  12*2 '3^1^  exP<-^2> 


-  exp 


7 


(A-3) 

which  corresponds  to  the  diagram  in  Fig.  4  and  is  one  of 
the  terms  contributing  to  (4-5) . 

For  the  discussion,  we  again  change  the  signs  of  the 
velocities  so  as  to  reverse  the  direction  of  motion  of 
the  vapor  molecules  in  Fig.  4.  Thus  ^>^2’ ^3  now  refer 
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to  the  initial  velocities,  \1 '  ,\2  ’  ,|3  ’  (|3  '  -£3)  to  the 
velocities  after  the  (12)  collision  and  "’^>3"  to 

the  velocities  after  the  (23)  collision. 

We  take  a  coordinate  system  in  which  molecule  2  is  at 
rest  at  the  origin  just  prior  to  the  (23)  collision  and 
the  Z-axis  is  in  the  direction  of  -f2'  (see  Fig.  8).  The 
velocities  of  the  molecules  in  this  coordinate  system  are 
at  various  times:  |i2,  “  ^±~i2'  ’  ^i'2'  and 

■+  4-  -)■ 

^ i ” 2 '  ~  ^  =  *>2,3).  The  center  of  the  droplet 

at  the  time  of  the  (23)  collision  (t  =  Ti2+T2'3^  is 

located  at  t  -  -p  -  f f12  -  I2'r2,3  =  -V12 -|2't2,3  for 

large  values  of  t12.  After  the  (23)  collision  molecule  2 

A  •> 

will  move  in  the  direction  of  C2,3  with  a  velocity  Z2"2 '  * 
In  order  for  molecule  2  to  catch  up  with  the  droplet,  the 
perihelion  vector  must  be  restricted  to  a  solid  angle 
proportional  to  r  -  |  +  ^2  T2T3f”  ’  Then  the  phase 

space  volume  behaves  as: 

/dT12dT2,3  /d^2'3  /dT12dT2’3l  ^lr12  +  ^2T2’3l 

(02)  (23)  (12)  (01) 


(A-4) 


which  diverges  as: 
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lim 

T'+'OS 

The  phase  space  volumes  associated  with  the  other  collision 

sequences  of  Pig.  3  also  require  the  occurrence  of  four 

successive  collisions  and  therefore  exhibits  the  same 

behavior  at  a  function  of  r^* 

In  addition  we  encountered  sequences  of  five  or  more 

collisions  among  three  molecules  and  the  droplet.  However, 

because  of  the  severe  restrictions  resulting  from  the 

required  occurrence  of  the  fifth  collision,  it  turns  out 

that  the  corresponding  integrands  fall  off  at  least  as 
_  2 

fast  as  t12  .  Thus  the  phase  space  volumes  of  these 

dynamically  more  complicated  collision  sequences  are 
finite. 


I 


dT 


12 


12 


=  lim  InT. 
T-*» 


(A- 5) 
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APPENDIX  B 

ANALYSIS  OF  HIGHER  ORDER  TERMS 


In  Section  IV  we  considered  an  expansion  for  the 
distribution  function: 

f  (?,e)  -  f(0)(p,I)  +  f(1)(p,1>  +  +  .  .  .  ,  (B-i) 

where  f^(p,£)  should  satisfy  equation  (4-2).  For  conveni¬ 
ence,  we  rewrite  this  equation  as: 


4-  S 2(p»012^1>^2)  +  S3^p,012’  J  *  (B-2) 


where 


-  ft0> 


<?,%)  + 
<p,i1)f<i)  (p,i2>- 


f(i)  (p,ipf(o)(?,y+ 


(B-2a) 


S2(p,B12,I1,e2)  (0)  (p,i3’)f(2)  (p,^)  +  f(2)  (p,^)f(0)  (p,lp+ 

-  f  (0)  (p^^f  (2)  Cp,%2)  -  f  {2)  (p/l^gCp,^)  , 

(B-2b) 
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with 


g Cp.b  -  f<0)(p,I)  -  1  (4), 


and 


S3(p,S12,l1,?2)  "  *a>  (p,ip  -  f(1)  (1)  (P,t2) 

This  equation  can  be  written  in  integral  form  as: 

j)  -  J<2)(?,!1)+a/'dT12  /d|2d312«:12e'aV<?l)Tl2 


’  [S2^“^1t12’  ®12’^1’^2^  +  S3^”^1T12,S12>^1’^2^  ’ 


(B-3) 


where 


u 

f (2)  (p,^)  =  a/dr12  Jd|2d6 
o" 


12*12® 


-av  (yr^ 


S1  (p"^lr12'  ®12*  *1'  *2*' 

(B-4) 


and  the  upper  limit  t  of  the  t^2- integration  is  again  given 
by  (3-4) . 

In  deriving  the  contribution  to  the  mass  flux  pro- 
2 

portional  to  a  lna  we  neglected  the  terms  S2  and  S3  in 

(B-2)  .  We  thus  identified  f ^  (p,^)  in  (4-9)  with  the 
~  (2)  ■* 

first  term  f  (p,^)  of  (B-3)  and  assumed  that  the  other 
term  in  (B-3)  would  not  modify  the  coefficient  of  this 
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logarithmic  term.  In  this  Appendix  we  shall  discuss  the 

evidence  for  this  assertion.  For  this  purpose  we  shall 

show  that  the  correction  term  to  f^  (p,^),  as- given  by 

(B-3) ,  leads  to  contributions  to  the  mass  flux  r  that  are 

2 

of  order  higher  than  a  lna.  The  mathematical  assumption 
implicit  in  this  procedure  will  be  discussed  at  the  end 
of  this  Appendix. 

The  correction  term  to  f ^  (p,|)  is  obtained  by 
iterating  (B-3)  once;  that  is,  we  approximate  f^(p,l) 

~  (n\  -+  •+ 

in  the  S2  term  by  f  (p,£).  This  then  gives  for  the 
2 

a  lna  contribution  to  the  mass  flux  T: 


-47rR2ran  a)  lim  —  /di|1£1*p  f^(p,£- 

a-0  a  lna 


-  +  <sf^ 


(B-5) 


Here  r 
*(2) 


(2) 


is  the  contribution  to  the  mass  flux  from 


(p,£)  and  is  precisely  the  coefficient  in  (4-12) .  The 
~  (2) 

correction  term  6?  may  be  written  as: 


6f(2)  =  6r{2)  +  5r^2)  , 


(B-6) 


where 
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*?»>  = 

f  (°)  (l’)f  (2)  (2’) 

1 

j(2)  (1»)f  (0)  (2’> 

+ 

- 

f(0)(l)  f(2)(2) 

L  J 

- 

f(2)(l)  g(2)  ]  , 

(B-6a) 

^  = 

f  (D  d')f  <!)  (2') 

- 

f (1)  f (2)  j 

,  (B-6b) 

and  we  have  again  used  the  shorthand  notation  introduced  in 
(4-13)  .  It  remains  our  task  to  show  that  and 

vanish  due  to  the  limit  cl+o  contained  in  these  terms. 

•wo)  —  (o) 

Considering  first  the  6r^  term,  we  decompose  f 
to  obtain: 


f(2>  -f}z,(5,li)  +  ?<«  Cp,l±), 


;(2) 


i(2) 


(B-7) 


where 


fp’cp.ii) 


a/dTiSL/d 


A 


:av^iTi3 


0 

0 


d^3dCi3Ki3e  *  “  Sl^-«iTi3>8i3^i'«3) 
for  e  I (p), 

for  ^  e  II  (p)  (B-7a) 


*ff  <?■«!> 


"  "av<^T13  ^  ->  ,  *  ■* 

dTi3  J  d^3d0  13*  i3e  S1  iTi3,C  i3  ’  ^  i  ’  ^  3^ 

0 


for  J,i  e  II  (p)  , 
for  e  I  (p) . 


(B-7b) 


and  T(p,^)  is  given  by  (3-5).  We  thus  have  in  analogy  to 
(4-15)  : 
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f^0)(r)f  j2)(2')  j  +  J2) 

f  J2)(r)fI(°)C2-)  ]+  [f^ 


f^}(nf  ]+  [fj(2) 

fi(0)(1)Jii)  (2)r  l?i2) 


<r)f{0)  (2*)  ]  +  | 

(l’)f  j2)(2')  ]+  [ 

(D  f  J  J^(2*)  J  -  [ 

(1) gj j (2)  j  , 


f^2)(l*)f  j0)(2’)]  + 

f  j0)  (l)f{2)  (2)]  + 
(B-8) 


where  gjj(2)  =  g{2)  gives  only  a  contribution  when 

f 2eI1  (|  P_ ^ iT  12 1  ^  we  now  sut)stitute  = 

and  f  ^  as  given  by  (2-14)  and  (3-6),  into 

(B-7a)  and  (B-7b) ,  and  then  in  turn  substitute  (B-7a)  and 
(B-7b)  into  (B-8) ,  we  find  that  6r^2^  is  given  by  a  sum  of 
integrals.  Each  of  these  integrals  is  associated  with  a 
collision  sequence  involving  five  or  more  collisions  among 
four  molecules  and  the  droplet.  The  leading  contribution  is 
determined  by  those  integrals  that  are  related  to  sequences 
of  five  successive  collisions;  these  sequences  are  presented 
in  Table  III. 

As  an  example  we  consider  the  collision  sequence 

(03) (34) (13) (12) (01)  which  is  schematically  represented 

~  (n} 

in  Fig.  9.  The  contribution  from  this  event  to 
is  given  by: 
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TABLE  III  Sequences  of  five  successive  colllelon*  among  four  eoleculea  end  the  dr  of  let 

associated  with  ST . 


Tern  in  (B.6) 

Collision  sequences 

■  ^£<«)(l')  f<*>(2')J 

(03) (34) (23) (12) (01) 

(04) (24) (23) (32) (01) 

(04) (34) (23) (12) (01) 

(04) (34) (23) (12) (01) 

(03)  (34)  (23)  (12)  (01) 

(02) (24) (23) (12) (01) 

(04) (53)(23) (12) (01) 

(02) (24) (23) (12) (01) 

(04)  (24)  (23X12)  (01) 

(02) (24) (23) (12) (01) 

(03)  (34)  (23)  (12)  (Cl) 

(04) (24) (23) (12) (01) 

(04) (24) (23) (12) (01) 

(04) (34) (23) (12) (01) 

(03) (34) (23) (12) (01) 

f<o)(2')] 

(03) (34) (13) (12) (01) 

(04)  (K)  (13)  (12)  (01) 

(04) (34)  (15) (12)(01) 

(04) (34) (13) (12) (01) 

(03)  (34)  (13)  (12)  (01) 

(01)  (14)  (13)  (12X01) 

(04) (34) (13)<12) (01) 

(01)(3l)(13)(12)<01) 

(04) (14) (13) (12>(01) 

(01) (14) (13) (12) (01) 

(03)  (34)  (13)  (12  X01) 

(04) (14) (15) (12) (01) 

(04)  (14)  (13)  (12X01) 

(04)  (34)  (13)  (12)  (01) 

(03)  (34)  (13)  (12)  (01) 

{f^U)  f(0(2)] 

(03)  (34)  (23)  02)  (01) 

(04) (24) (23) (12) (01) 

(04)(34) (23)  (12)  (01) 

(04)  (34)  (23)  (15)  COl) 

(03)  (34)  (23)  (15)  (01) 

(02)  (24)  (23)  (U)  (01) 

(04)  (34)  (13)  (12)  (01) 

(02 >(24) (23) (12) (01) 

(04)  (24)  (23)  (12)  (01) 

(02>(i4) (23)  (12) (01) 

(03) (34) (23)  (15) (01) 

(04)  (24X23)  (12X01) 

(04>(24)  (23  )  02)  (01) 

(04)  (34)  (23)  ("5)  (01) 

(03)  (34)  (23)  (12X01) 

[fCo)(l<)  f(2)<2')l 

l  I  II  -■ 

(02) (24){23) (12) (01) 

(01)  (l4)  (13)  (12X01) 

[f<o)(l)  f<^(2)] 

(02) (24) (23) (12) (01) 
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(2)  2 

6r  (03)  (34)  (13)  (12)  (01)  "  ~47r  R  mn~cu  lim 


a 


“  ro  a-0  lna 


00  00  00 


/dT  12  7dTl,3  Jdr 3’4  l/d^ld^2d^12d^3d°l’3d^4d°3,‘ 


0  0  0 


A  n  2 

K12Kl,3K3T4^l'^  exP{-av<Ti2 +Ti,3)}  exp{-(^ +£  ”  )} 


^  exp  {  -  ( V + 


•f2 


»2 


-exp 


f«s’a  +  ■ 


(B-9) 


where  we  have  replaced  v(£^)  and  v  (£  j)  by  the  constant  v 
since  this  change  will  not  affect  the  result.  The  meaning 
of  the  velocities  ^  and  f /’  is  explained  in  Fig.  9. 

The  symbols  t. 0.,  and  jc..  refer  to  the  collision  that 
produces  the  velocities  and  £  ^ •  In  order  to  study  the 

behavior  in  the  limit  cr+0,  we  need  to  estimate  the  phase 
space  for  large  values  of  and  T^.g.  Using  the  same 

arguments  as  in  Appendix  A  it  is  evident  that  Bg,4 
is  restricted  to  a  solid  angle  proportional  to 
l^lr12  +  ^1t1'3  +  ^3T3’4l  *  After  carrying  out  the 
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integration  over  Tg,^,  the  integrand  will  be  proportional 
to  +  ti^rsf-1  Thus  the  asymptotic  behavior  of  the 

integrand  is  determined  by  an  integral  of  the  type: 


(B-10) 


Using  the  coordinate  transform  x  =  avTi2  Rnt*  y 
we  see  that: 


I 


I?!*  +  fj'y| 


J 


avr^,3  , 


(B-ll) 


which  vanishes  in  the  limit  a-*  0.  Similar  arguments  can 

be  used  to  show  that  the  contribution  to  (B-6)  from  all 

sequences  in  Table  III  vanishes  in  the  limit  a *  0. 

If  we  next  consider  the  term  in  (B-6b) ,  we  find 

that  it  may  also  be  decomposed  into  a  sum  of  integrals, 

each  of  which  may  be  associated  with  a  collision  sequence. 

These  collision  sequences,  however,  involve  six  or  more 

collisions  among  four  molecules  and  the  droplet.  The  phase 

space  volumes  associated  with  these  sequences  are  thus 

smaller  than  those  associated  with  the  sequences  in 
~  (2) 

Table  III  and  will  likewise  be  zero. 

In  the  main  body  of  this  report  we  extracted  from 
2 

equation  (4-2)  an  a  lna  contribution  to  the  mass  flux  r 

~  (o) 

and  determined  its  coefficient  r  .  In  this  Appendix, 
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we  have  demonstrated  that  the  leading  contributions  from 

the  terras  neglected  are  of  higher  order  in  a  and  do  not 
~  (O') 

modify  r  .  A  completely  rigorous  proof  of  this  logarith¬ 
mic  term,  however,  should  consider  all  of  the  higher  order 

terms  and  show  that  their  sum  is  also  of  higher  order  than 
2 

a  lna.  For  that  purpose  we  would  not  only  need  the  solu¬ 
tion  of  (4-2)  to  all  orders  in  a,  but  also  the  sum  of  the 
higher  order  terms  in  (4-1) .  In  this  respect,  the  situation 
is  again  analogous  to  the  density  dependence  of  the  trans¬ 
port  properties  of  gases.  For  a  three  dimensional  gas  a 
re summation  of  the  most  divergent  terms  encountered  in  a 

formal  power  series  expansion  in  the  density  n  leads  to  the 

2  ~ 
n  In  n  term  in  (5-2)  with  a  finite  coefficient  r)2  [16, 

20,22] .  Again  for  a  rigorous  proof  one  would  have  to  show 

that  a  resummation  over  all  the  less  divergent  terras  would 

not  change  this  result.  Such  a  complete  proof  is  not  yet 

available  for  the  density  dependence  of  the  transport 

properties,  nor  for  the  Knudsen  number  dependence  of  the 

mass  flux  considered  in  this  report. 
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